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In Pedersen et al. pp we proposed a method to utilize a temporally dependent superlattice potential 
to mediate spin-selective transport, and thereby probe long and short range magnetic correlations 
in optical lattices. Specifically this can be used for detecting antiferromagnetic ordering in repulsive 
fermionic optical lattice systems, but more generally it can serve as a means of directly probing 
correlations among the atoms by measuring the mean value of an observable, the number of double 
occupied sites. Here, we provide a detailed investigation of the physical processes which limit the 
effectiveness of this “conveyer belt method”. Furthermore we propose a simple ways to improve the 
procedure, resulting in an essentially perfect (error-free) probing of the magnetic correlations. These 
results shows that suitably constructed superlattices constitute a promising way of manipulating 
atoms of different spin species as well as probing their interactions. 

PACS numbers: 37.10.Jk, 73.21.Cd, 74.25.-q, 75.50.Ee 


I. INTRODUCTION 

Gold atoms in optical lattices are powerful “quan¬ 
tum simulators” for strongly correlated many-body sys¬ 
tems in periodic potentials. The experimental study of 
these systems has resulted in a range of breakthrough 
results in the last two decades, including the observation 
of the Mott-superfluid transition ma, fermionic pair¬ 
ing Dirac points [B], and various topological band 
structures [IH]. Fermionic atoms in optical lattices 
can be accurately described by the single Fermi Hub¬ 
bard model [M [U, and a major goal of current re¬ 
search in this area is to explore this model at low en¬ 
tropies where it, e.g., is expected to exhibit antiferro¬ 
magnetic order. Moreover, there is strong numerical evi¬ 
dence that the Fermi Hubbard model captures the essen¬ 
tial physics leading to high-Tc superconductivity [IM7]. 
Using Bragg spectroscopy [18] , one recently observed the 
onset of antiferromagnetic correlations reaching temper¬ 
atures as low as 1.4 times the critical temperature for an¬ 
tiferromagnetic ordering m, and such correlations were 
also observed along one dimension in an anisotropic lat¬ 
tice HO] • Several alternative methods have been proposed 
for observing magnetic correlations, including time-of- 
flight pTM2d| . and light-polarization experiments [5ill25j . 
A common problem for most of these probes [2(11124] is 
that they are only indirectly probing, e.g., the antifer¬ 
romagnetic ordering. Since they inherently rely on dy¬ 
namics in momentum space [lainHis] or on averaging 
over large areas |25j , they are not very sensitive to direct 
spatial correlations of the spins such as the antiferromag¬ 
netic ordering of neighboring atoms. 

In [T] , we demonstrated how a moving superlattice po¬ 
tential can be used like an atomic conveyor belt to se¬ 
lectively transport one atomic spin component a desired 
number of lattice sites. In this way, one directly can 
measure the spatial dependence of the spin-spin corre¬ 


lation function, and in particular detect both the long 
range order as well as short range correlations which may 
mark the onset of long range correlations. Importantly, in 
this proposal the desired correlation function is directly 
mapped into an observable quantity, the number of dou¬ 
ble occupied sites. This means that there is no need to 
average over several shots of the experiment since correla¬ 
tions can be directly measured in a single shot. This will 
be a particularly important property if the details of the 
experiment is fluctuating from shot to shot so that not 
every instant of the experiment show the desired order¬ 
ing. In this article we investigate the scheme of Ref. [1] 
in further detail, improving and optimizing the different 
steps. 

A detailed description of the superlattice conveyer belt 
(SCB) probing scheme of Ref. [1] is presented in Sec. 

An essential step in the procedure is the velocity with 
which the superlattice is moved through certain avoided 
crossings to ensure the correct adiabatic or diabatic be¬ 
havior in order to obtain the desired selective hopping of 
atoms. We analyze this in Sec, jll B] by mapping the sys¬ 
tem to a Landau-Zener tunneling between nearly degen¬ 
erate levels on neighboring sites. Based on this mapping 
we identify suitable parameter regimes where the proce¬ 
dure can have a good performance. Then in Sec. m we 
consider unwanted hopping between lattice sites. Again 
we reduce this to a simple model by identifying it with 
decoherence caused by the width of the Bloch bands. 
Based on the simple model we find suitable regimes of 
superlattice amplitudes where we predict the scheme to 
have a good performance. In Sec. m we consider an¬ 
other source of erroneous hopping, due the presence of 
resonances of higher Bloch bands as the superlattice is 
moved. These hoppings are a major source of error for 
the original procedure of Ref. [1] , but we show that they 
can be almost eliminated by ramping down the superlat¬ 
tice amplitude in between the desired hoppings. Finally, 
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in Sec. Ill we investigate the validity of our conclusions 
from the simplihed models with a numerical calculation. 
These simulations show that the modified procedure can 
measure the spatial dependence of the spin-spin correla¬ 
tion function essentially without error both for short and 
long distances. 


II. THE SUPERLATTICE CONVEYER BELT 
A. Ideal procedure 

In the following, we assume that the dynamics in the 
transverse y and z directions is frozen by a tight trans¬ 
verse lattice potential, so that the problem is effectively 
one-dimensional involving only the cc-direction. We con¬ 
sider atoms trapped in an optical lattice described by a 
potential of the form 

Va{x) = (I) 

where A is the wavelength and A denotes the strength 
of the potential. With this potential it is convenient to 
introduce the lattice recoil energy, Er = h'^/2m\^ as 
a natural energy scale while time is measured in units 
of [TD]. As a specihc example, we consider 

atoms in two different spin (hyperfine) states, but the 
method may be generalized to atoms with more spin 
states. At time t = 0 the atoms are prepared in the 
many-body state we wish to probe. How we arrive at 
this is not important for the analysis below and it may, 
e.g., be an equilibrium state or the result of a dynami¬ 
cal evolution of a strongly interacting system. From the 
time we apply the SCB probe, however, we assume that 
the interactions among the atoms are sufficiently weak 
that we can ignore them. This can be the case if the in¬ 
teractions are sufficiently weak that they can be ignored 
on the time scale of the SCB, which can be much faster 
than the preceding dynamics, or if the interactions are 
turned off before the SCB using for instance a Feshbach 
resonance pS] . 

The first step of the probing procedure is a spin- 
selective Raman pulse which excites one of the spin 
states to the hrst excited vibrational level, while leav¬ 
ing the other spin in the vibrational ground state m- 
A time-dependent superlattice potential with twice the 
wave length of the form 

VsixA) = B{t) sin^ {7r[a; - XB{t)]/\} (2) 

is then applied. Here B denotes the strength of the su¬ 
perlattice potential and XB{t) is a time dependent dis¬ 
placement of the superlattice potential, which we use to 
create the conveyor belt below. This potential needs to 
be turned on sufficiently fast that atoms do not have 
time to tunnel to neighboring sites, but at the same time 
sufficiently slowly that they adiabatically follow the vi¬ 
brational levels at each site. 


The key idea in our probing procedure is to make the 
time-dependent potential act as a “conveyor belt” for 
atoms populating the first excited vibrational level. To 
do this we use the total potential Va + VB{t) to create 
an array of almost isolated double well systems. This 
happens at xb ^ A/4, where the maxima/minima of the 
superlattice potential VB{t) are close to the maxima of 
the original lattice Va, as illustrated in Fig. (a). We 
shall henceforth refer to this position as the degeneracy 
point. At this position atoms can not easily tunnel be¬ 
tween lattice sites of Va, which have a maximum of Vb 
in between them. On the contrary tunneling is not sup¬ 
pressed between sites which have a minimum of Vb in 
between them. Thus effectively the lattice potentials cre¬ 
ates an array of coupled double wells. In these double 
wells the atoms in the vibrationally excited states move 
between degenerate levels located in the two wells, while 
the atoms in the vibrational ground state will remain al¬ 
most stationary due to the higher tunneling barrier. To 
turn this into a conveyor belt we need to interchange the 
populations of the first excited states of the two wells. 
We do this by moving the superlattice potential across 
the degeneracy point. In Fig. (a) the superlattice Vb 
is displaced slightly to the left of the degeneracy point. 
This means that the right well is slightly higher in en¬ 
ergy. Moving the superlattice to the right exchanges the 
roles so that the left well is now higher in energy. If the 
superlattice is moved slowly enough the populations will 
follow adiabatically and thus be interchanged. Due to the 
difference in the tunneling strength one can tune the su¬ 
perlattice velocity vb = dxs/dt such that only the atoms 
in the first excited vibrational level follow adiabatically, 
whereas atoms in the ground state will remain station¬ 
ary as illustrated in Fig. (a). This corresponds to the 
atoms in the first excited states following the instanta¬ 
neous Bloch level spectrum adiabatically, as opposed to 
atoms in the ground state, which cross the level spectrum 
diabatically. This is illustrated in Fig. [^around t = lOOt/j 
(see below for how we obtain this plot). 

The above procedure moves atoms in the first ex¬ 
cited vibrational level by a single position. To make the 
conveyor belt move the atoms over larger distances, we 
quickly displace [Fig.[^(b)] the superlattice potential Vb 
to the next degeneracy point where a new array of double 
wells are created, but now with each well coupled in the 
opposite direction [Fig. (c)]. Here it is essential that 
the superlattice is moved quickly in-between the degen¬ 
eracy points, since we need to ensure that atoms remain 
at the same positions during this period. In a Bloch band 
picture this corresponds to the diabatic level crossings il¬ 
lustrated around e.g. t « and t « 2Qt)EB, in Fig. [2| 

and is discussed in detail in Sec. |HD[ Repeated crossings 
of degeneracy points then move the atoms in the first 
excited state down the lattice. The number of crossings 
then directly translates into the number of translations 
of the atoms. Note that atoms on odd (even) sites move 
to the left (right) during this procedure. 

At the end of the procedure, the superlattice is turned 
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off and the number of doubly occupied sites is measured. 
This can for instance be done by transferring atoms 
from the excited vibrational level to the ground state 
and merging the doubly occupied states into molecules 
through a Feshbach resonance or photo association 
m [28] followed by a measurement of the number of 
molecules. The fraction of molecules to atoms then di¬ 
rectly measures the long-range spin-correlations in the 
system as discussed in Ref. [T|. Specifically in the ideal 
limit the molecular fraction, Pmoiijn)-, after m applica¬ 
tions of the conveyor belt, can be related to the spin 
correlation function through, 

Pmol = ll^-{sU^st). (3) 

Here (sf+mSf) is the spin-spin correlation function be¬ 
tween spins, located m sites apart. 

In total the SCB is created by a number of different 
sequences with different parameters. The time evolution 
of the various parameters entering into the evolution is 
shown in Fig. (a) for the particular example analyzed 
in Ref. [T|. In the following, we analyze in detail the 
different steps of the SCB procedure in order to minimize 
the main sources of error. We furthermore improve the 
method by allowing the amplitude B of the superlattice 
potential to vary with time. We show that this may 
significantly reduce the error of the protocol so that it 
can operate with almost no error. 

B. Optimizing the superlattice velocity 

The crucial step in the SCB procedure involves the 
neighbor well degeneracy shown in Fig. [^a) and (c) and 
Fig.i where all atoms in the first excited vibrational 
state move to their neighboring site, while all atoms in 
the ground state remain at their initial site. If however, 
the superlattice moves too slowly close to the degeneracy 
point, atoms in the ground state will also tunnel to their 
neighbor sites. If, on the other hand, the superlattice is 
moved too fast, atoms in the first excited state do not 
have sufficient time to tunnel to the neighboring site. It 
follows that there is an optimum superlattice velocity, 
where the degeneracy point is traversed neither too slow 
nor too fast, in order to ensure that most of the atoms 
in the first excited state move to the neighbor site, while 
most atoms in the ground state stay put. 

We will now use Landau-Zener theory for adiabatic 
passage in two level systems to estimate the optimum 
velocity. Our system, however, does not consists of two 
level systems and we need to extract effective param¬ 
eters to use in the Landau-Zener theory from the de¬ 
scription of the full system. Regardless of the position 
of the superlattice, the instantaneous potential is always 
periodic with a period A set by the periodicity of the 
superlattice. This means that we can use Bloch’s theo¬ 
rem which specifies that the eigenstates with energy E\. 
of a one-dimensional periodic system are Bloch waves, 
(a;|<I>],) = v}f.{x), parametrized by their wave number. 



odd even odd even odd even 



FIG. 1. Snapshots of the superlattice “conveyer belt” pro¬ 
cedure. (a) Superlattice potential Vb (red line) placed near 
the degeneracy point where neighboring wells are degenerate. 
The combined potential creates distinct double well struc¬ 
tures, e.g., counting from the left, minima number one and 
two of Va (blue line) are connected, whereas the tunneling 
from minima two to three is suppressed by Vb- Here the su¬ 
perlattice velocity is tuned such that atoms in the first vibra- 
tionally excited state (red circles) cross adiabatically moving 
along with the conveyor belt (thick arrows), whereas atoms 
in the ground state (blue circles) remain stationary. Atoms 
in the odd sub-lattice move in the opposite direction of atoms 
starting in the even sub-lattice. The dashed arrows show po¬ 
tential Landau-Zener errors, while wiggly arrows show delo¬ 
calization errors due to the presence of multiple double wells, 
(b) After the crossing, the superlattice is moved quickly to the 
next degeneracy point. Note that atoms starting in the even 
sub-lattice are only weakly confined allowing them to delocal¬ 
ize easily as indicated by the stronger wiggly lines, (c) The 
next degeneracy point, which transports atoms in the excited 
states to the next neighboring site. Here a double well system 
is formed by, e.g., minima number two and three of Va- (d) 
Modified superlattice procedure with a decreased superlattice 
strength during the displacement. In the modified SCB this 
step replaces part (b). 


k, and their vibrational level, I- We therefore calculate 
the Bloch level spectrum numerically for various lattice- 
super lattice amplitudes {A,B) and displacements xsit)- 
This is done by numerical diagonalization of the Hamilto¬ 
nian discretised on a set of 256 real space positions within 
the unit cell. The result of this calculation is shown in 
Fig. i (plotted as a function of time rather than dis- 
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FIG. 2. The original superlattice procedure with A = 40i?fl 
and B = 20Eji proposed in Ref. [I], (a) Superlattice-lattice 
parameters as a function of time, when the superlattice con¬ 
veyor belt moves atoms two sites away from their starting 
point. Note that v is small but non-zero in the intervals 
around t — lOOtij and t = 2501/?, whereas it is exactly zero 
around t = 0 and t = 330ti? (b) The instantaneous Bloch 
eigenenergies for all Bloch wavevectors fc as a function of time. 
The Bloch eigenenergies for each wavevector almost coincide, 
so that the Bloch bands look like lines. The dashed arrows 
indicate the possible Landau-Zener errors. Zoom of the eigen¬ 
energies near the avoided crossing for (c) the ground state and 
(d) the first excited state. The superlattice velocity is chosen 
such that the ground states cross diabatically, while the first 
excited states cross adiabatically. 


placement). Note that while the figure seems to consist 
of discrete lines, each line is in fact a whole band of en¬ 
ergies corresponding the various fc-values, but the width 
of each band is too narrow to be visible. 

Close to the degeneracy point, the plots in Fig. (c) 
and (d) show that the system is well described by two 
level systems corresponding to tunneling between one vi¬ 
brational level from each well. Furthermore they have 
nearly linear asymptotics far away from the crossing, i.e. 
E 2 — El = (3x as required by the Landau-Zener descrip¬ 


tion. At the avoided crossing, we approximate the system 
by two states, of energies Ei and E 2 , tunnel coupled by a 
matrix element set by the gap A between the levels at the 
degeneracy point. We extract the energy gap A and the 
slope P from the energy spectrum for each Bloch wave- 
vector. The probability of an atom moving diabatically 
from one instantaneous eigenstate to the other through 
the avoided crossing is then given by the Landau-Zener 
transition amplitude EiEni 

Plz = (4) 

Near the degeneracy point, the level spacing is pro¬ 
portional to the displacement of the superlattice and 
we thus have a = vbP- A perfectly working SCB will 
have the ground state eigen-energies crossing diabati¬ 
cally, Phzigs) ~ 1, while the first excited state energies 
cross adiabatically, Plz (ex) « 0. The transition proba¬ 
bility can then be calculated from Eq. @ and the total 
success probability when crossing the degeneracy point 
is defined as the average over all the Bloch wave vector 

^crossing = ^ ^ '\ PLzigs) — Plz{s-x)]. (5) 

k 

From this the optimum superlattice velocity is found by 
maximizing Pcrossing- 

In Fig. we plot the resulting optimal error probabil¬ 
ity, 1 —Pcrossing for lattice amplitudes {A, B) in the range 
[4Ai?, 40Aj?]. Note that this estimate requires that the 
Landau-Zener description is a sufficiently good approx¬ 
imation of the dynamics. We therefore have to exclude 
an area in the lower left corner defined heuristically by 
A > 5/4, where we find that the curves do not have near 
linear asymptotics and is thus not well described by the 
Landau-Zener formula. This is, however, anyway not the 
regime we are interested in, since the scheme does not 
have a good performance for these parameters. On the 
other hand in the regime where the performance is good, 
we find the Landau-Zener approximation to be justified. 
From the figure we identify a region in (A, B) space where 
the error is less than one percent at the optimal veloc¬ 
ity. This requires a large A, and B has to be larger than 
a minimal value which is nearly independent of A but 
smaller than an upper limit set by A. 

Below we find a more stringent requirement for the 
lower limit for B and we shall therefore not discuss it 
further. The requirements on the value of A has a sim¬ 
ple physical explanation: A large A is required in order 
to have a large energy separation between the ground 
state and the first excited state, and therefore a large 
difference in the tunneling rates of the two levels. On 
the other hand a strong potential Vb pushes the minima 
of the double well potential closer to the barrier in the 
middle and thus effectively lowers the tunneling barrier. 
In fact for B > 4A the bump in the middle of the dou¬ 
ble well disappears. Hence in order for the system to 
be meaningfully described as a double well, there is an 
upper limit on B which grows with increasing A. 




































5 


40 

35 

30 


iu‘25 

U) 

1 20 
CQ 

15 

10 

5 

5 




FIG. 3. Logarithmic color plot of the error probability of the 
Landau-Zener procedure as obtained from Eqs. Q-([^. The 
white dashed line encloses the region with an error less than 
one percent. The hatched area marks parameters where the 
dynamics cannot be described by Landau-Zener physics. 


C. Bloch level delocalization close to the neighbor 
well degeneracy 

The SCB procedure is based on the fact that at the 
degeneracy points, the lattice can be approximated by 
an array of isolated double wells. This approximation 
is, however, not exact and atoms can still hop between 
different double wells which leads to delocalization errors 
as indicated by the green arrows in Fig. This delocal¬ 
ization is most significant for the first excited vibrational 
mode, and can be used to construct a simpler but also 
less controlled probe of ordering in optical lattices [3T]. 
For the SCB probe, on the other hand, the delocalization 
is a source of error and in this section, we shall investigate 
the effect of it. 

In order to ease our analysis, we label the states of the 
Va lattice by a superlattice site index n and a compound 
index I = {m,r]), with m referring the vibrational level 
and 77 = e/o to the even or odd sub-lattice. Note that this 
notion of Wannier states being localized in the original 
lattice Va is strictly speaking only meaningful away from 
the degeneracy point, but this suffices for our present 
purpose. 

The Wannier wave functions localized around the n’th 
superlattice site can be constructed from the instanta¬ 
neous Bloch waves, 

<(ai, to) = ^^ to). (6) 

Loosely speaking, the delocalization error is due to the 
the spread in energy of the Bloch band which means that 
the different wave numbers have different phase evolu¬ 


tion. This difference in phase evolution is just the math¬ 
ematical description of the atoms tunneling to different 
lattice sites. While the width of the Bloch band may 
seem negligible near the degeneracy point in Fig. one 
must remember that it should be compared to the slow 
pace set by the optimum Landau-Zener velocity. 

We define the localization probability as the over¬ 
lap between the initial first excited state (ffi), and 

the adiabatically time-evolved state, which is 

transferred from the neighbor well, 77 , during the crossing 
of the degeneracy point. We will show that this amounts 
to an integral over the energies within the Bloch band. 
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(7) 

( 8 ) 


To see this, we first write down the adiabatically time- 
evolved state, 


Due to the spatial symmetry of the initial and final state 
the instantaneous Bloch eigenstates are related in the 
following way, Combining 

those two facts allow us to turn Eq. Q into Eq. ([^ with 
only a bit of algebra. 

From the localization probability, we define the local¬ 
ization error, 1 — Piocai- As seen from Eq. ^ it depends 
on the time available for particles to tunnel to a different 
site. To extract this time, we need to consider how wide a 
region around the degeneracy point we need to traverse. 
The superlattice velocity is slowed down in a range of su¬ 
perlattice positions xb G [A/4 — 7i;/2,A/4-|-w/2], where 
suitable widths w are determined heuristically: Assum¬ 
ing that the potential minima do not change significantly 
from the position of the minima of the potential Va the 
derivative of the energy is /3 ^ 2 ttBIX. Hence by setting 
w = 2AX/B the energy separation at the endpoints is 
E 2 — El ^ 2 ttA which is much larger than the splitting 
A at the avoided crossing. Combined with the optimum 
superlattice velocity obtained in the previous section, we 
use this to calculate the transversal time ti — and 
then the delocalization error for a variety of the lattice- 
superlattice amplitudes {A,B). 

The results are shown in Fig.|^ As seen, the delocaliza¬ 
tion error decreases with increasing superlattice potential 
strength, B because a stronger superlattice reduces the 
coupling to the neighboring sites. The delocalization er¬ 
ror is less than 1% above the solid line. We also indicate 
by a hashed line the region where the Landau-Zener er¬ 
ror is less than 1% as also shown in Fig. In this way, 
the solid and hashed lines define a region in {A, B) space 
which is well suited for the SCB procedure, since both 
the delocalization error and the Landau-Zener error is 
less than 1%. 
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FIG. 5. Eigenenergies of the original superlattice procedure 
with A = 40E_r and B = 30Er. Arrows show the desired path 
taken by the ground state and first excited state starting from 
the even lattice sites. The filled grey circles mark unwanted 
resonances between the ground state, first excited state, and 
higher excited vibrational states. 


FIG. 4. The delocalization - as defined by Eq. 0 - displayed 
on a color log scale as a function of lattice amplitude A and 
superlattice amplitude B. The delocalization is less than one 
percent above the solid line. The dashed line is the same as 
that in Fig.j^and shows where the Landau-Zener crossing has 
an error probability of less than one percent. The hatched 
area marks parameters where the dynamics cannot be well 
described by Landau-Zener physics. 


D. Decreasing errors between neighbor well 
degeneracies — the modified SCB 

The two previous sections analyzed the possible errors 
associated with the traversal of the degeneracy point, and 
determined a range of values for the lattice parameters 
{A,B), where the SCB procedure works reasonably well. 
Now, we discuss the dynamics when the superlattice is 
in-between degeneracy points, as illustrated in Fig. 0b). 
We will first show that certain complications arise during 
this step in the original SCP procedure and then discuss 
how the procedure can be modified to yield an almost 
perfectly working conveyor belt by turning down the su¬ 
perlattice potential, when it is in-between the degeneracy 
points. 

As seen from Fig.|^ the superlattice amplitude B must 
be comparable to the lattice potential, A in order to 
avoid the delocalization of the first excited state during 
the traversal of a neighbor well degeneracy. A large su¬ 
perlattice potential, however, causes errors when moving 
between consecutive neighbor well degeneracies. During 
this step, we want all the atoms to remain at the same lat¬ 
tice position but delocalization errors as well as unwanted 
resonances with higher excited states in the neighbor well 
can diffuse the atomic wave functions. Examples of such 
unwanted resonances are shown in Fig. where the un¬ 
desired resonances are marked with grey circles. These 
errors are particularly large during the displacement if a 
strong superlattice potential is used since the superlattice 


pushes one of the states towards being almost unbound 
as shown in the right hand side of Fig. (b). This can 
directly be seen by comparing Fig. with Fig.[^ The lat¬ 
ter has a higher value of B and hence more and stronger 
avoided crossings. 

The unwanted effects discussed above can be decreased 
by speeding up the superlattice potential, ensuring that 
the wave function has no time to delocalize and that any 
avoided crossings, associated with the unwanted reso¬ 
nances, are traversed diabatically [T]. The superlattice 
velocity cannot, however, be increased indefinitely, since 
this induces uncontrolled excitations to higher Bloch 
bands if the states cannot follow the bands. Hence one 
must find optimal crossing velocities which balances in 
between both sources of errors. For parameters such as 
the ones used in Fig. the avoided crossing are so strong 
that we have not been able to find a suitable velocity 
where the SCB has a good performance. For the orig¬ 
inal proposal of Ref. [1] we have therefore only been 
able to identify a small region of parameter space around 
{A,B) = (40Efl,20i?fl) where the procedure works well. 
Even for these parameters, however, the error is consid¬ 
erable. For atoms in the first excited state occupying the 
sub-lattice site which is lifted in energy by the superlat¬ 
tice [the red atom on the right hand side of Fig. Iil(b)l 
the error grows to approximately 60% after 7-8 repeti¬ 
tions [T]. 

Fortunately, the source of error discussed above can be 
almost eliminated by using a modified superlattice “con¬ 
veyor belt” (MSCB) method, which only uses a strong 
superlattice potential in the vicinity of each degeneracy 
point. By ramping down the superlattice in between the 
degeneracy points there is strongly reduced tunneling of 
the atoms to unwanted sites. The ramped down super¬ 
lattice potential is shown in Fig.j^d), which is the MSCB 
replacement of Fig. [^b) of the original SCB. 

As a tangible example, Fig. 0a) shows the time- 
evolution of the superlattice position, xb, the su- 
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FIG. 6. The modified superlattice conveyor belt procedure 
with neighbor well degeneracies traversed at ^ = 407?^ and 
B = SOEr. (a) The lattice-superlattice parameters as a func¬ 
tion of time, (b) The instantaneous Bloch energies as a func¬ 
tion of time. Note how the intermediary ramp down of the 
superlattice potential, B, removes unwanted energy crossings 
present in the original procedure of Fig. 


perlattice velocity, vb, and the superlattice intensity 
for the modified procedure when neighbor well degen¬ 
eracies are crossed at lattice-superlattice amplitudes 
{A, B) = {A0Eii,30EB). Both the superlattice ampli¬ 
tude and the superlattice velocity is changed smoothly 
in order to decrease unwanted excitations. Figure ib) 
shows the instantaneous eigen-energies of the superlat¬ 
tice Bloch waves as a function of time for the MSCB. 
In contrast to the original SCB procedure of Fig. the 
unwanted crossings of the Bloch bands have now been 
removed by the dynamical ramp down of the superlat¬ 
tice amplitude. This eliminates the main sources of er¬ 
rors in-between degeneracy points, and the accuracy of 
the MSCB method is now mostly limited by the physics 
of the adiabatic/diabatic evolution near the degeneracy 
point. 


III. NUMERICAL CALCULATIONS 

In order to quantitatively test our predictions con¬ 
cerning the MSCB, we numerically calculate the single¬ 
particle dynamics of the system by solving the time- 


dependent Schrodinger equation. As explained above, 
we assume that the atomic wave functions remain in 
the vibrational ground state in the perpendicular (y, z)- 
directions. This assumption is valid if the total time 
taken to perform the superlattice procedure is much 
smaller than the wave function dephasing time along 
these directions, which can be realized experimentally 
by using a tight lattice potential in the (y, z)-directions. 
Working in superlattice units (i.e. measuring x in terms 
of A, energies in terms of Er and times in units of tu), 
the one dimensional Schrodinger equation describing the 
system is then 

( 10 ) 


We discretize the Schrodinger equation on a real-space 
grid with 128 points and solve the dynamics by expand¬ 
ing a a basis consisting of the 16 lowest lying eigenstates 
of the underlying lattice Va- The MSCB procedure is 
pieced together from basic steps, 1... z... m, defined over 
a range of superlattice positions xb € [xBii),XBii + 1)]. 
The steps fall into three categories: the initial and final 
superlattice ramps, the dynamics close to the degeneracy 
point, and the dynamics in-between degeneracy points. 
The dynamics is calculated as evolution matrices U for 
each of these basic steps of the MSCB procedure, as dis¬ 
cussed in Appendix This allows for an efficient cal¬ 
culation of the total evolution simply by multiplying the 
evolution matrices for each individual step together. In 
the original SCB, the superlattice velocity and superlat¬ 
tice intensity were kept constant within each step, but as 
illustrated in Fig.[^ they were allowed to change abruptly 
between consecutive steps. In contrast, the MSCB only 
allows the superlattice parameters to vary smoothly in 
time. Specifically, we chose the superlattice parameters 
to be third order polynomials in time with vertical tan¬ 
gents at the end points of each basic step. As explained 
in Appendixj^this amounts to a total of 9 parameters in¬ 
cluding the optimum superlattice velocity, which can be 
changed in order to find the optimal MSCB procedure. 

We choose the lattice-superlattice amplitudes at the 
degeneracy point, and use the optimum superlattice ve¬ 
locity calculated from Eqs. as a reasonable start¬ 

ing point. The remaining parameters are chosen so as to 
maximize the total localization probability of the atoms 
in the first excited states, estimated as 


pi — 
local 


1 

N ^ 


( 11 ) 


In order to measure the quality of the MSCB we define 
the error probability as the amount of a wave function 
which does not end up in the desired place. The ground 
state wave functions, / = 1, should remain put, while 
the first excited state wave functions, I = 2, should be 
moved a number of lattices sites, m, corresponding to the 
number of traversed degeneracy points, m. 


































IV. CONCLUSION AND DISCUSSION 



\i 

a 


u 

M 


1 

i_ y _M 


1 _M 


k a 

Excited state 

.... even sites 



d * 



U 

Ground state 

.... even sites 
odd sites 

. . . . i.i . . 


-5-4-3-2-1012345 
X (units of X) 


FIG. 7. The ground state and excited state Wannier func¬ 
tions of even and odd lattice states after 1-8 repetitions of 
the MSCB procedure where the degeneracy point is traversed 
with lattice-superlattice amplitudes {A,B) = (40i5ij, SOFlij). 


Fixing the lattice-superlattice amplitudes {A, B) at the 
degeneracy point, we perform multiple simulations of the 
MSCB each time varying one or more of the remaining 
parameters. Because of its vastness we cannot map out 
the entire parameter space, and we are satisfied when we 
find a local minimum. 

In this way, we find that the MSCB procedure 
works best for degeneracy point lattice strengths 
{A,B) = {A0E]i,30Eji), where it is essentially without 
error. This is illustrated in Fig. where we plot the 
wave functions after up to eight lattice translations with 
{A,B) — {A0E[i,30Eji). We see that the vibrational 
ground state wave function is essentially unperturbed, 
whereas the first vibrational states are translated up to 
eight lattice positions as desired. 

The error probability of each of the four relevant ini¬ 
tial states obtained for various amplitudes of the lat¬ 
tice and the superlattice is plotted in Fig. The 
excellent performance of the MSCB procedure with 
{A,B) = {A0Efi,?>0Efi) is evident. In addition, many of 
the other scenarios perform as well as the best example 
of the original SCB procedure described in [T]. As can be 
seen, the simulated errors diminish with increasing val¬ 
ues of A or B, which is in qualitative agreement with the 
error estimates of Fig. and Fig. Only the extreme 
case {A, B) = {AOEn, AQEr) abhors from this picture, be¬ 
cause the large superlattice potential makes it impossible 
to separate the different levels when moving between de¬ 
generacy points. 


We have carefully analyzed the different steps in the su¬ 
perlattice conveyor belt and identified the limiting factors 
and potential sources of errors. Based on this, we have 
been able to optimize and improve the performance by 
dynamically tuning the superlattice potential during the 
procedure. The main features of our hndings obtained 
by simple model calculations, were confirmed by a direct 
simulation of the system dynamics for various lattice- 
superlattice parameters. Based on the optimization, we 
found an almost ideal performance of the conveyor belt 
for lattice parameters around {A,B) = {A0Eji,30Eii). 
As a result the probe can be used to perform an es¬ 
sentially error free measurement of magnetic correla¬ 
tions in only a single shot of the experiment by mea¬ 
suring the fraction of doubly occupied sites. As opposed 
to experiments based, e.g., on measuring noise correla¬ 
tions [2I1E21E5] where an average has to be performed 
over multiple shots, this is a major advantage of this 
proposal. We thus believe that the procedure can be a 
valuable tool for future optical lattice experiments. 

From an experimental perspective one of the main 
challenges of the proposed method is that it requires a 
precise phase lock between the lattice laser potential and 
the superlattice laser potential. Under the ideal condi¬ 
tions (A, B) = {iOEfi, 30Eii) we have assumed that the 
cross over region in the adiabatic passage has a width of 
w ~ 0.009A. Hence the relative position of the two lat¬ 
tices needs to be controlled with great precision. Possibly 
this requirement can be reduced by working with higher 
excited vibrational states or by extending the region of 
the slow dynamics. The latter will, however, come at the 
prize of a longer operation time of the procedure. An¬ 
other challenge for the proposed method arises from the 
fact that in most experimental systems the lattice height 
is not homogeneous but varies across the sample. Hence 
the lattice parameters cannot be kept at ideal working 
conditions for all positions in the lattice. Here it is a 
major advantage of our modified scheme that it is much 
more robust to variations than the original scheme which 
only worked in a small region in parameter space. In the 
end the procedure is based on Landau-Zener transitions 
which are not very sensitive to the precise value of the 
tunneling. As long as the parameters do not show too 
strong variations, we therefore believe that it should be 
possible to apply the modihed procedure even for systems 
which are not completely homogeneous. 
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FIG. 8. The error probability, given as the wave-function weight left outside the desired final state, for the MSCB procedure 
for various lattice amplitudes A, and superlattice amplitudes B. Each plot shows the error-probability as a function of lattice 
translations for atoms in the ground vibrational state and in the first excited state. Note that the MSCB procedure distinguishes 
between atoms on even and odd lattices sites as determined by the initial position of the superlattice. The lattice parameters 
with the smallest error are {A, B) = {AOEr, SOEr), for which the central crossings can be optimized almost perfectly, and the 
disturbances in-between these crossings can be almost eliminated. The parameterizations of each of these MSCB scenarios are 
explicitly stated in Appendix [A| 
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Appendix A: Parametrization of the MSCB 

Here we describe the detailed parameterization of the 
MSCB, which was used to obtain the results in Figs. 
and The parameters and basic steps are illustrated 
in Fig.|^ We parameterize the modified SCB procedure 
in terms of five different control points. The first point 
controls the initial step I with a slow ramp up of the su¬ 
perlattice, where the intensity is increased to some small 


value Bi over a period ti. The initial ramp is followed 
by step II, where the superlattice is further ramped up 
while being simultaneously moved towards the next de¬ 
generacy point. The superlattice position is advanced at 
a high velocity, V 2 , until it reaches the second control 
point, after which the velocity is smoothly ramped down 
in order to the match the optimal velocity, V 3 , upon the 
arrival at the third control point. Specifically we model 
the velocity ramp as a third order polynomial in time 
with horizontal tangents at the initial and the final con¬ 
trol points. At the same time the superlattice intensity is 
increased further. This is done smoothly by following a 
third order polynomial in time with a horizontal tangent 
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b) Parametrization 

I. Initial ramp and t, 

II. Intermediate ramp [d^VJ^, and p = X3-X2 

III. Neighbor well deg. and w = x^-X3 

FIG. 9. a) The parametrization of the modified SCB proce¬ 
dure. The hve control points for the three different stages of 
the MSCB: the initial ramp, the intermediate ramp and the 
neighbor well degeneracy. Note that the procednre is sym¬ 
metric around the degeneracy points. Circles mark the five 
control points and dashed lines indicate the horizontal tan¬ 
gents at the end points of the third order polynomial curves, 
b) Summary of the parametrization of the MSCB for each of 
the three stages. 

at the initial control point. 

At the second control point, we choose a superlattice 
intensity B2 and a slope [dB/dt]2- The ramp continues 
smoothly in time from the second to the third control 
point and ends at a superlattice intensity, B3 - again 
with a horizontal tangent. We chose step II to happen 
over a range of superlattice positions, p = X3 — X2 (see 
Fig.§, during a total time, T = 2p/(ua — ui). 

At the next step III, the degeneracy point is traversed 
at a constant superlattice intensity, B^, and velocity, U 3 . 
After traversing the degeneracy point we have another 
step ir, where we decrease the superlattice intensity and 
increase the velocity. This is done in a completely sim¬ 
ilar fashion as the ramp before the degeneracy point in 
step II by constructing the time reverse of that procedure 
with the same parameter values. This ensures that the 
procedure is symmetric with respect to the degeneracy 
points at superlattice positions xb = A/4 + nXj2, where 
n is an integer. 

As described in the main text, we choose fixed values of 
the lattice-superlattice amplitudes (A, B) at the degener¬ 
acy point. We then estimate the optimal Landau-Zener 
velocity by optimizing Eq. The remaining parameters 
are initially chosen such as to minimize the overall delo¬ 
calization error as estimated by Eq. From this initial 
point in the impossibly large parameter space, we per¬ 
form multiple simulations where we vary each parameter 


slightly in order to find a local minimum of the total error 
probability. 

In Figurej^we showed the error probability of the best 
found parameterization of the SCB for different lattice- 
superlattice amplitudes (A, B) at the degeneracy point. 
We supply the parameterizations for each (A, B) pair in 
the table below: 


(A,B3) 

[Er] 

Vbi 

[Er] 

ti 

[tR] 

V2 

[A/ifi] 

Vb2 

[Er] 

diVeb 

[ER/tll] 

P 

[A] 

W 

[A] 

10^«3 

[A/tfl] 

(20, 20) 

2.0 

2.0 

0.5 

6.00 

35.3 

0.1 

0.03 

150 

(20, 30) 

2.0 

0.4 

0.5 

6.00 

54.3 

0.1 

0.10 

150 

(30, 20) 

0.5 

1.0 

0.1 

6.00 

9.55 

0.1 

0.01 

2.0 

(30, 30) 

6.0 

6.0 

0.15 

6.47 

3.92 

0.2 

0.03 

4.5 

(30, 40) 

1.0 

1.0 

0.05 

6.00 

5.85 

0.08 

0.035 

40.0 

(40, 20) 

1.0 

10.0 

0.005 

9.00 

0.95 

0.08 

0.004 

0.12 

(40, 30) 

6.0 

30.0 

0.015 

8.90 

1.00 

0.16 

0.009 

0.35 

(40, 40) 

8.0 

8.0 

0.1 

8.95 

1.00 

0.2 

0.012 

0.30 


Appendix B: Evolution matrices 

In this appendix we describe the calculation of evo¬ 
lution matrices, which forms the basis of our numerical 
simulation. Consider the matrix solutions U{t,0) to the 
ordinary partial differential equation 

if = A{t)U, (Bl) 

where A{t) = A(t + T) is a periodic time-dependent ma¬ 
trix with a period T. The equation x = A[t)x with 
the boundary condition a;(0) = afg is solved directly 
by x{t) = U{t, 0)xq if we impose the initial condition 
17(0,0) = 1. Applying the periodicity of A we see that 
the U matrix can be decomposed into 

U{t + mT,0) = U{t,0)U{T,0)"^, m £ N. (B2) 

In our simulation we solve the evolution within each 
Bloch wave sector separately, from 

$^(t) = (B3) 

where k is the Bloch wave number, and p. is a level index. 

Because the Bloch bands do not mix, we find the general 
evolution matrix elements between different levels, p and 
u, for each wave number fc, 

= (B4) 

In some common basis, |$^), we have 

and they are naturally concatenated by matrix multipli¬ 
cation in order to yield the transformation matrix for the 
compound evolution, 

Uln{t2,to) = Ulm{t2,tl)Umn{tlAo)- (B5) 
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Each step of the procedure is then solved independently, 
multiplied together and extended to multiple periods 
through the use of Eq. |B2| The final evolution can be 
transformed to a real-space position basis by rewriting 
U^n ™ 9' basis of Wannier functions, 


a at site i at time t. The workings of the SCB is then 
described by transition matrix elements, defined by 

= (Cl) 

j 


(B6) 

VNk ^ 

such that the real-space evolution matrix becomes. 


Write a = a(to). The probability that the site i is doubly 
occupied after the sweep so a molecule can be formed is 
fmoi = where Ui = djd^. In the Mott 

phase we have that 


^ k 


(B7) 


Choosing the common basis, to coincide with 

the initial (and hnal) eigenstates of the system, we can 
directly identify the probability that an atom in some 
initial Wannier state end up another Wannier state after 
the time evolution, 

= (bs) 

The success probabilities can be calculated as the norm 
square of the appropriate matrix elements. Starting in 
the ground-state n = 0 at a position Rq = 8, this atom 
should remain put and the success-probability is given 
by Psuccess = |t7oo(^)P' When instead starting in the ex¬ 
cited state /To = 1, the atom should move along the lat¬ 
tice, and the success-probability after an even number of 
repetitions, m, is given by Psuccess = 


Appendix C: Correlation functions 

The correlation functions measured by the SCB can be 
determined in the following way. At an initial time to we 
begin the application of the SCB, which ends at a later 
time, ti- The operator dl^^{t) creates a particle with spin 


Pmoi = E 

± E . (C2) 

Introduce the spin operators = l/2-|-s^, = 1/2 —s^, 

s’*' = and s~ = Then 

Pmoi = E - (s|sf)) 

31 

± E {spT). (C3) 

Here we assumed (s^) = 0. Knowing the transition ma¬ 
trices one can easily calculate the double occupancy if 
knowing (s|sf) and {s'j sf) from e.g. Quantum Monte 
Carlo simulations. 

When the spin down species is pinned we have 
^— ^ij 5 S-nd 

Pmoi = E - (s|sf)). (C4) 

3 

If the spin up species is transferred perfectly m sites along 
the lattice, U^ij = and 

Pmoi = 1/4 — (C5) 
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